Perturbation of infinite networks of resistors 
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The resistance between arbitrary nodes of infinite networks of resistors is studied when the network 
is perturbed by removing one bond from the perfect lattice. A connection is made between the 
resistance and the lattice Green's function of the perturbed network. Solving Dyson's equation 
the Green's function and the resistance of the perturbed lattice are expressed in terms of those of 
the perfect lattice. Numerical results are presented for a square lattice. Our method of the lattice 
Green's function in studying resistor networks can also be applied in the field of random walks as well 
as electrical and mechanical breakdown phenomena in insulators, thin films and modern ceramics. 
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I. INTRODUCTION 

The jcjalculation of the resistance between two arbitrary nodes of infinite networks of resistors is a well studied 
subjectE □. Recently, it has been demonstrated how the krttice Green's functions, as an alternative method to using 
the principle of the superposition of current distributionsclQ, can be applied to this problemU. The interested readers 
will also find several useful references there. 

The analytical behavior pitithe lattice Green's function has been extensively studied in condensed matter physics 
over the past three decadesElQ; impurities are often modeled by a simple forms of perturbations (see Eq. (13) below). 
Thus, the resistor network problem can be successfully tackled by utilizing the strong connection between the two 
fields, in particular through the use of the lattice Green's functionEl. Our approach may also serve as a didactic 
introduction of the Green's function method used in solid state physics as well as field theory. The perturbation series 
that show up in these fields are, in general, infinite and need to undergo a tedious renormalization procedure. In 
contrast, the perturbative solutions of the Dyson equation can be summed up in closed form in our case, due to the 
special form ( 13) of the perturbation. Studying the analytical behavior of the closed- form Green's function has also 
proven pedagogically useful in introductory courses on solid state physics and field theory. 

The conductivity of lattices with randomly distributed defects has also been thoroughly investigated in the past 
25 years-within the framework of effective medium theorieaS'B and by using the position-space renormalization-group 
methodtB. Below, we shall use the Green's function approach to study a special case, i.e. to determine the resistance 
for a so-called perturbed lattice that is obtained by removing one bond from the perfect lattice (corresponding to 
"simple perturbation" mentioned in the previous paragraph). Although this system is much less complex than those 
studied previously, it has the great advantage of being analytically treatable and leading to closed-form expressions 
for the resistances. 

As an example (see Fig. n^), consider an infinite square lattice whose edges represent identical resistances R. Re- 
moving one edge (bond) from this perfect lattice results in a perturbed lattice. 




FIG. 1. Perturbation of an infinite square lattice by removing one edge between sites Tig and Vjg. Each edge represents a 
resistance R. The resistance R{i,j) is measured between arbitrary lattice points r^ and Tj. 

Then, one can ask the resistance between two arbitrary nodes. It is simple to fpjd the resistance between the two 
ends of the missing bond. On one hand it is well known (see e.g. Aitchison's papeiO) that for a perfect square lattice 
the resistance between adjacent nodes is R/2. On the other hand this resistance equals the parallel resultant of R 
and the resistance we wish to find. Thus, it is obvious that the resistance between the two ends of the missing bond 
in the perturbed lattice is R. 

Finding the resistance between two arbitrary nodes in the_perturbed lattice is much more difficult. In this paper we 
present a further application of the Green's function methodO to answer this question. Our treatment of the problem is 
based on the idea used by Kirkpatrickl for studying the transport in inhomogeneous conductors. The Dirac notation 
of the vectors is used throughout this paper. Note that the derivation and our formulas presented below can be 
straightforwardly expressed by the usual matrix-element formalism. However, the Dirac notation seems to provide 
the most powerful and elegant formulation of the problem. The readers who arer-iinfaniiliar with this formalism will 
find sufficient background material in several well-known textbooks, e.g. SchwabliS. 

According to Ohm's and Kirchhoff's laws, for a given current distribution, the potentials on the lattice sites are 
given by a Poisson-like equation involving a so-called lattice Laplacian operator. This operator can be decomposed 
into two parts: one associated with the perfect lattice and the other corresponding to the perturbation. A relation 
between the resistance and the Green's function for the perturbed lattice is derived. The problem thus ceduces to the 
calculation of the Green's function for the perturbed lattice which actually satisfies the Dyson equationQ. As we shall 
see, the Dyson equation can be solved exactly in the present problem. Then, an explicit formula can be derived for 
the resistance in the perturbed lattice in terms of the resistances between different nodes in the perfect lattice. Our 
Green's function method proved to be a highUt-effective technique and can successfully tackle the present problem 
even in cases when the superposition principlecrQ faces hardly (or even in-)surmountable difficulties. 

II. GREEN'S FUNCTION AND RESISTANCE FOR PERFECT LATTICE 

To establish the method and the notations that we shall use throughout this paper, the case of a perfect lattice is 
reviewed first. Consider a perfect d-dimensional infinite lattice made up of identical resistances R. All lattice points 
are specified by position vectors r given in the form 

r = ?iai+/2a2H h/^ad, (1) 

where ai, a2, • • • , a^ are independent primitive translation vectors, and Zi, I2, ■ ■ ■ ,ld range through all integer values. 
We denote the potential at site r^ and the current entering at lattice point r^ by V{ri) and /(r^), respectively. 
According to Ohm's and Kirchhoff's laws, we may write 

J2[V{r^)-V{r, + n)]=RI{r,), (2) 

n 

where the n are the vectors from site r^ to its nearest neighbors. One can form two vectors, V and / from the 
potentials and the currents at sites r.;: 

i 

I = Y.\')h, (4) 

i 

where Vi = V{Yi) and li — I{vi). Here it is assumed that |j), associated with the site r^, forms a complete orthonormal 
set, i.e. {i\k) = 6ik and J2i l*)(*l — 1- Using the vectors V and /, Eq. (||) can be rewritten as 

E(z%-A,,)(j-|l^ = i?(z|/, (5) 

3 

where z is the number of neighbors (coordination number) of each lattice site (e.g. z = 2d for a d-dimensional 
hypercubic lattice), and A^i is unity if the sites r^ and r; are nearest neighbors and zero otherwise. The summation 
is taken over all lattice sites. Premultiplying both sides of the above equation by \i) and summing over i, we have 



L^V = -RI, (6) 

where Lq is the so-cahed lattice LaplacianH: 

Lo = El*)(^^^--^'5,,)(j|. (7) 

For further progress it is useful to define the lattice Green's function of Lq- Similarly to the common definition 
used in the literatureEl, the lattice Green's function is defined by 

ioGo = -1. (8) 

The properties of this lattice Green's function and its explicit form in coordinate representation are given in Appendix 
^. The solution of the Poisson-like Eq. ^ can then be given in a simple form: 

V = -RL^^I = RGqI. (9) 

To measure the resistance between two arbitrary sites we assume that a current /p enters at site r^ and a current 
— /q exits at site Vj, while the currents are zero at all other lattice points. The current distribution may thus be 
written as 

Irn = la {Srni - S^j) for aU m. (10) 

The potentials at different sites can be determined by inserting the above-given current distribution into Eq. M), and 
one finds 

Vk = {k\V = R{k\GoI = Rj2iMGoH Im = RIo[Goik, i) - Go(fc, j)] , (11) 

m 

where Go{n, m) = {n\Go\m) is the matrix element of the operator G in the basis |n). The resistance between sites r^ 
and Tj is then 

Roii,j) = -^^V^ = 2R[Go{i,t) - Go{i,j)] , (12) 

-'0 

where we have made use of the symmetry properties of the Green's function Go(i, j) given in Appendix H. The same 
expression for the resistance in terms of the lattice Green's function was derived in Ref.El. 

III. GREEN'S FUNCTION AND RESISTANCE FOR PERTURBED LATTICE. DYSON'S EQUATION 

The current contribution Sit at site r^ due to the bond (ioja) is given by 

= (*l(Ko)-bo))((*o|-Oo|)^=(*|ii^, 
where the operator ii is of a so-called 'dyadic' form: 

ii = (No)-bo))((*o|-Oo|). (13) 

Here we used the fact that (5„,„ — {n\m). Removing this bond from the perfect lattice the current li at site r^ is given 
by 

{~LoV).^-R5h = Rh. (14) 

Thus, Ohm's and Kirchhoff 's laws for the perturbed lattice can be written as 

LV = -RI, where (15) 

L = Lo + Li. (16) 



Equation (|15|) is similar to that found in the case of a perfect lattice. However, the operator L is now a sum of 
Lq associated with the perfect lattice and a 'perturbation' given by Li. Note that the same decomposition of the 
Laplacian operator with the perturbation Li given by (|l^) was used by Kirkpatrick. 

Similarly to the perfect lattice, the Green's function G for the perturbed lattice is defined by 

LG=-l. (17) 

The resistance for the perturbed lattice can be obtained in terms of the Green's function much in the same way as 
for a perfect lattice. Measuring the resistance between sites r^ and Yj, we assume that a current /o enters at site 
Yi and a current — /q exits at Vj. The currents are zero at all other lattice points. Thus, the current distribution 
is /„ = /o {5rm - Smj) for all m. From Eqs. © and (0) we have Vk = {k\V = R(k\GI = R(k\GJ2„^ l™> Im = 
R'^j^{k\G\m) /„ — i?X)m ^i^j '^) ^rn- Inserting the current distribution /„ one finds 

Vk=RIo[G{k,t)-Gik,j)]. (18) 

Therefore, the resistance between sites r^ and Tj reads 

Rihj) = -^^V^ = R[G{i,^) - G{i,j) + G{j,j) - G{j,i)]. (19) 

Jo 

Note that unlike for a perfect lattice, G{i, i) ^ G{j,j) since translational symmetry is broken in the perturbed lattice. 
However, as we shall see, G{i,j) = G{j,i). The problem of finding the resistance reduces to the calculation of the 
Green's function for the perturbed lattice. p. 

Perturbation theory for the Green's function is worked out in the literature (see e.g. Economou's bookQ). Using 
Eqs. (P), ( p^ and ( p^ ) yields (— G,^^ + Li) G — —1. Premultiplying both sides of this equation by Gq one obtains 
the so-called Dyson's equation: 

G = Go + GoLiG. (20) 

This is an equation for G in terms of Gq (which is assumed to be known), and the perturbation Li. In general, the 
solution of the Dyson equation can be found by iteration resulting in an infinite series: 

G = Go + GqLiGq + GqLiGqLiGq + • • • . (21) 



However, if Li has a special form - such as that given in Eq. (13) - the Dyson equation can be solved exactly. Now 
the perturbation Li is equal to the dyadic product of the vector \io) — \jo) with itself. In this case one may apply the 
identity 

valid for arbitrary vectors \x) and \y) whose dimensions are the same as that of operator A, assuming that the inverse 
of operator A exists and the denominator 1 + (j/|A~^|a::) ^ 0. This identity is readily proved by postmultiplying Eq. 
( P^ ) with A + \x){y\. Using the above identity with A = Lq, \x) = |io) — \jo) and {y\ ~ (io| — (jol one obtains for the 
Green's function: 

^_ (r ,r \-^-n , Go (No) - Ijo)) (fal - (jol) ^o .„.,^ 

1 -2[Gq(iq,io) -Go(io,jo)J 

Below we shall see that the denominator in Eq. ( p3| ) is never equal to zero when d > 1. The matrix elements of the 
operator G can then be expressed with the matrix elements of Go as 

nf ■\ /-ir-iA nr -^ , [Go(»,»o) - Go(^, jo)] [Go(^o, j) - Go(jo, j)] 

G{i,j)^{i\G\j)^Go{i,3) + . — n . (24) 

l-2[Go(io,«o) -Go(«o,Jo)J 

There is an alternative way to calculate G. This method is analogous to that presented in Chapter 6. of Economou's 
bookQ for the calculation of the Green's function. There the system was condensed matter modeled by a periodic 
tight-binding Hamiltonian with one impurity given by Dirac delta potential. 



One can insert Li given by (13) into the series in (21). Then the matrix elements of G result in infinite geometric 



series, which can be summed up exactly leading to the same results as above. 



Note that this geometric series as well as the perturbative solution of the Dyson equation given by Eq. (|21| ) is not 
necessarily convergent, especially in the case of dyadic perturbations with |(i/|^~^|a;)| > 1. However, Eq. ( p2[) is still 
vaHd even for |(y|A~-'^|a;)| > 1. 

It is clear from Go{i,j) — Go{j,i) that G{i,j) is also a symmetric matrix, i.e. G{i,j) — G{j,i). ^From Eqs. (|l9| ) 
and (EJ) the resistance between r^ and Tj can be obtained in terms of the matrix elements of Gq: 

:^=G(*,z) + G(j,j)~2G(*,j) 

= 2[Go{t,t)-Go{i,j)] + -^ — — — — — ^ '—. (25) 

1 - 2[Go(«o,»o) - Go(«o,Jo)J 

Using Eq. (112) this formula can be rewritten in terms of the resistances in a perfect lattice: 

R{h3) = Ro{hj) + -F^ ^—. — -Ti ■ (26) 

4:[R- Ro(io,]o)\ 

This is our final result for the resistance between two arbitrary nodes r^ and Tj of the perturbed lattice in which the 
bond (io jo) is removed. 

It is efl,sy to calculate the resistance between sites r^g and r^g for a d-dimensional hypercubic lattice. For symmetry 
reasonscl the resistance between r^Q and Tj^ in a perfect lattice is i?o(«0: jo) = R/d if d > 1. Then, using Eq. (|2q), 
the resistance between the two ends of the removed bond is i?(«o, jo) = R/{d — 1). For a square lattice {d — 2) the 
resistance is R as mentioned in the introduction. It also follows that the denominator in Eq. (Eq) is never equal to 
zero, and this is necessarily true for each previous equation that contains this term. 

Note that the explicit form of the lattice Laplacian Lq defined in Eq. (M) was not used in the derivation of Eq. (pQ). 
Thus, our final result for the resistance in the perturbed lattice is valid for any lattice structure in which each unit 
cell has only one lattice site. This is the case, for example, for triangular lattices. _ 

Finally, it is worth mentioning some possible applications of our method and results. In Doyle and SncU's bookEI 
the connection between random walks and electric networks is presented. Thus, our result can also be extended to the 
random walk problem in a perturbed lattice. Furthermore, the highly effective method of the lattice Green's function 
in studying resistor networks also provides a pedagogically useful way to familiarize students with the notions of the 
Dyson equation and the Green's function. 

When more than one bond is missing from the perfect lattice, our method (outlined above) can be iterated. It can 
be shown that our final result ( pq ) is still valid provided the values of Ra of the perfect lattice are replaced by the 
values of the resistances obtained from the previous step of the iteration. Thus, lattices with more complex defects 
can be studied analytically with our method. An example is the so-called crack-type defects arising in-several fields - 
such as electrical and mechanical breakdown phenomena in insulators, thin films and modern ceramical3. The current 
distribution around cracks can be easily calculated from the Green's function. This probleiiuwas studied numerically 
as well as analytically in the continuum limit (for a recent paper see e.g. Boksiner and Leathllj and references therein). 
Application of our method to investigate this problem is in progress. 

IV. NUMERICAL RESULTS 

Below we shall present some numerical results for square lattices. Recurrence relations for the resistances between 
arbitrary nodes were derived for perfect square lattices in Ref.El. These recurrence formulas provide a very simple and 
effective tool for the calculation of the resistance. For convenience, we invoke them: 

Am 2m — 1 „ . 

Rnlm + l,m + 1) ~ Rn{m,m) Rnim — 1, m — 1), 

^ ' ^ 2m+l ^ ' ^ 2m+l ^ '' 

Ro{m + 1, m) = 2RQ{m,,m,) — Ro{m, m — 1), 
Ro{m -1-1,0)= 4:Ro{m, 0) - Ro{m - 1, 0) - 2Ro{m, 1), 
Ro{m + l,p) = ARQ{m,p) - RQ{m-l,p) - Ro{m,p+l)- Ro{m,p-l) if < p < m, (27) 

where the indices m and n of Ro{m, n) (not to be confused with i and j above) stand for the relative coordinates of the 
two sites on the square lattice. Since the exact values of i?o(l, 0) = R/2 and i?o(l, 1) = 2R/'k are knownS (obviously 



i?o(0, 0) = 0), one can calculate the resistance exactly for arbitrary nodes by using the above recurrence formulas. It 
is interesting to note that using these formulas one has, in units of R 

6400711399252571342562758751832284129928 1 _ 
°^ ' ' ~ 1089380862964257455695840764614254743075 n ~ ' 

which can be found quickly (within a few second) using the program MAPLE. In Fig. Hthe resistance Ro{i,j) between 
the node i = (0,0) and j — {jx,jy) is plotted as a function of j^ and jy for a perfect square lattice. 




FIG. 2. The resistance Ro{i,j) between nodes i — (0,0) and j — {jx,jy) as a function oi jx and jy for a perfect square 
lattice. 

One can see fronr the figure that by increasing the distance between the nodes i and j the rotational symmetry of 
the resistance is more and more recovered. Indeed, in Ref.U it was shown that for large values oi jx or/and jy, the 
asymptotic form of the resistance is 

R (, r-, ^ ln8 



^0- Jj.) = -V^ \hl + Jy' + 7 + — ) , (28) 



where 7 = 0.5772 ... is the Euler-Mascheroni constant. 

On the perturbed square lattice the resistance can be calculated from Eq. (pq). As an example we show results 
when site r^ is fixed and Vj is moved along the line of the removed bond. The resistance is measured between nodes 
Vi and Tj, where r.^ = (0, 0) and Vj — {j^, 0). The two ends of the removed bond are r^^ = (0, 0) and rj^ = (1, 0). In 
Fig. H the resistances for the perfect and the perturbed lattices are plotted as functions of j^ ■ 
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FIG. 3. In an infinite square lattice the origin of the coordinate system is at site r^, i.e. r^ — (0, 0). The ends of the removed 
bond are at Vig — (0, 0) and rj„ = (1, 0). The resistances on the perturbed and the perfect square lattices are measured between 
sites Ti and r-,, where r-, — {jx,0)- They are plotted as functions of jx (filled circles and empty squares, respectively). The 
solid and dashed lines are just for guiding the eyes. 

One can see that the resistance is always larger in the perturbed lattice than in the perfect lattice. This intuitively 
obvious fact follows from the positivity of the second term in Eq. (|2g). It can also be seen from the figure that 
increasing the distance between nodes r^ and r^ the resistance tends to that of the perfect lattice. The resistance is 
not symmetric under the transformation jx — > ~jx because translational symmetry is broken in the perturbed lattice. 
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FIG. 4. The same as Fig. H except that the ends of the removed bond are at ri„ — (1, 0) and rj^ — (2, 0). 

The configuration is the same in Fig. ^ as in Fig. y except that r^,-, = r,; + (1,0). Similar features are observed in 
the two cases. However, the resistance on the perturbed lattice tends to that of the perfect lattice more rapidly in 
Fig. U than in Fig. || as jx increases. 

Finally, it worth mentioning that a similar analysis of the resistance can be carried out for simple cubic lattices. 
Recently, Glasser et al.li3 have expressed the exact values for the cubic lattice Green's functions in terms of gamma 
functions. Using their results one can calculate the resistance for perfect simple cubic lattices and for perturbed 
lattices. 
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APPENDIX A: EXPLICIT FORM OF THE LATTICE GREEN'S FUNCTIONS FOR PERFECT LATTICE 

STRUCTURES 

In this section we shall present the symmetry properties of the lattice Green's function Go, and derive its explicit 
form in coordinate representation for perfect lattice structures. 
From Eq. (Q) we have Lo(i, j) = («|^o|j) = Lo{j,i). Therefore, 

Go{t,j) = {i\Go\j)^Go{j,i). (Al) 

Thus, Go{i,j) is a symmetric matrix. The perfect lattice is translationally invariant, so Go{i,j) depends only on the 
relative coordinates of sites r^ and r j . Thus we have 

Go{i,i) = Go{j,3). (A2) 

Let us turn to the derivation of the coordinate representation of Gq. Consider a finite part of the lattice constructed 
by repeating the unit cell (say Ni times) along the unit vectors a^, where i = 1,2, ■ ■ ■ ,d. Define the Fourier transform 
of the basis \i) by 

|k)-;^E^''^""'l*)- (A3) 



The lattice sites r^ are specified by Eq. (M, and the vectors k belonging to the first Brillouin zone (BZ) are given by 

where mi, TO2, • • • nid are integers such that — 7Vi/2 < rrii < Ni/2 for i = 1,2,- ■ ■ ,d, and the hj are the reciprocal 
lattice vectors defined by a^bj = 27rSij, i,j = 1,2, ■■■,(!. The total number of lattice sites is A^ = N1N2 ■ ■ ■ N^. 
Periodic boundary conditions are assumed at the boundary of the finite lattice, and in the final results the limit 
A?^ ^ 00 is taken. It is also assumed that Ni is an even integer, which will be irrelevant in the limit N ^ 00. The 
mathematical description of crystal lattices and the concept of the Brillouin zone can be found in many books on 
solid state physicals. To distinguish the basis vectors |i) from their Fcutfier transforms |k), parentheses will be used 
instead of angled brackets. Using the so-called lattice summation rulestj 

^g»(k-k')r, ^^^^^,^ (A5) 

i 
^g,:k(r,-r,)^^^^^.^ (A6) 

k 

it can be shown that the basis |k) is a complete orthonormal set, i.e. (k|k') = (5kk' and ^j^ |k)(k| = 1 provided the 
basis \i) is complete and orthonormal. It is also clear that 

(z|k) = -^e''^'-'. (A7) 



The Laplacian operator Lq for a perfect lattice is a diagonal matrix in the basis |k) since 

(k|io|k') = J2(Mi) (A,, -z6,,) (j|k') = ^J2 ^"'*'"" (^^^- -'^^^) ^"'"^''^ = -^(1^) -^kk', (A8) 



A^ 



where 



e(k) = z-^e'''", (A9) 

n 

and vectors n point from some lattice point into its the nearest neighbors. 

In coordinate representation the matrix element of the Green's function Gq can be written as 

Go(^,J) = (z|Gob-> -5](z|k)(k|Go|k')(k'|^) ^ 1 ^ (k|Go|k')e^(''-'^'"-^^ (AlO) 

kk' kk' 

In the third step the completeness of the basis |k) has been used. Since Go = —Lq, one finds 

1 ,-^ e'''(r>-r3) 
kGBZ ^ ' 



Note that the toroidal and cylindrical geometry cases discussed by JenglJ can be straightforwardly treated by direct 
summatiDn over k. For infinite lattices (in the limit N — > oo) the discrete k summation can be substituted by an 
integratj: 

V ^ ^ / d% (A12) 

kGBZ (27^) -^keBZ 

where w^nis the volume of the unit cell of the lattice. Thus, the lattice Green's function in coordinate representation 
becomeaJ 

JkeBZ (27r) £lKj 



One can see that Go(«,j) depends only on the relative coordinates of sites r^ and r^, as stated in Eq. (A2). For 



example, z = 2d for the d-dimensional hypercubic lattice, and one finds e(k) = 2 ( d — ^^^^^ coskai j . For body 

centered cubit and face centered cubic lattice structures see e,g. the introduction of the lattice Green's function by 
Katsura et alE, and for triangular and hexagonal lattices Ref.cl. 
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